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Tlii« study oviilualcMl the* upplicatiou of LnndsAt multinpcctiral scanner 
(MSS) imagery for assessing the trophic status of the major iakes and 
reservoirs in the Tennessee Valley, 

Ground truth water quality data collected by the Environmental 
Protection Agens^y during the 1973 National Eutrophication Survey (NES) 
at 35 reservoirs in the greater Tennessee Valley region were subjected to 
cluster and principal component data analyses to develop a trophic state 
index for the reservoirs. Each reservoir’s trophic state index was 
defined by its relative position on a first principal component axis. 

Water quality characteristics selected as trophic indicators used in the 
cluster un.ilysis and principal component analysis were clilorophyll a, 
conductivity, total phosphorus, total organic nitrogen, the inverse”of 
the Eecchi disc depth, and the yield of an algal assay procedure. 

Lindsat MSS data from four different dates were extracted from com- 
puter tapes using a semi-automated digital data handling and analysis 
system at the University of Wisconsin - Madison. Reservoirs were extracted 
from the surrounding land matrix by using a Band 7 density level slice of 
3; and descriptive statistics to tnclude mean, variance, and ratio between 
bands for each of the four bands were calculated. 

Significant correlations (>0,80) between the MSS statistics and many 
trophic indicators were identified. Regression models were developed to 
predict reservoir eutropivieation using MSS statistics as the independent 
variables and the tropbic state inticx, developed from ground truth, as 
the dependent variable in each frame. Regrosiiion models were also devel- 
oped to predict Secchi disc depth, conducUivity, and total phosphorus. 

The models gave significant estimates of each reservoir's trophic state 
as defined by its trophic state index and explained in all four tandsat 
I i’v(me,s at least 85 pet cent (R^) of tl»e variability in the data. Each 
l«nidsat I fame bad its own unique models which were not practically 
applicable on other dates, 

To illnstralc the spatial viiriaLluns within reservoirs as well as 
the relative variations between reservoirs, a table-look-up elliptical 
classificatioiv was used in conjunction with each reservoir’s trophic 
state index to classify each reservoir on a pixel -by-pixel basis and 
produce color-coded thematic representations. 

Although the need for ground truth information on water quality 
places a restriction on the use of Landsai MSS data for the prediction 
of trophic state, such data still has value in regions where there arc 
many lakes or reservoirs wUhln a s ingle MSS frame. Under such cir- 
cumstances the collection of ground truth from a small number of ’'bench 
mark” lakes for the development of regression models would result in 
considerable cost savings. 
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INTRODUCTION 


The fcience of remote lensingi thit ii the collection of informetion 
about an object from a distance, is a rapidly developing water quality 
monitoring tool that is expected to complement in situ sampling and 
analysis. Remote sensing has rapidly progressed from a purely qualitative 
water quality monitoring tool to one which permits quantitative prediction 
of certain water quality characteristics over large areas based on limited 
ground observations. Some rather obvious benefits of remote sensing arc 
the ability to present a synoptic pictorial representation of an extensive 
area as opposed to a specific location, to suirvey vast areas in a very 
short time, and to observe at a later time information that was not fully 
realized or being sought at the time the data were collected. 

Host past Landsat eutrophication studies (Boland, 1976; Scarpace 
et al. , 1976) have focused on natural lakes. In this study, reservoirs 
were investigated. The relatively short hydraulic retention times of 
most reservoirs ensure that responses to changes in waste loads will be 
detected in a matter of years rather than the decades required for most' 
naturally formed lakes. Host nutrient loading trophic state models have 
been developed using data from natural lakes in the northeastern United 
States and Europe (Vallenweider, 1968; Rast and Lee, 1978). Reservoirs 
have relatively high flushing rates and respond differently to nutrient 
inputs than do natural lakes with relatively lower flushing rates. 

In this study we have investigated the application of Landsat multi- 
spectral scanner imagery to assess water quality conditions in the 
Tennessee Valley region, with particular emphasis on the determination of 
the trophic status of the major reservoirs in the Valley. The purpose of 
this vroject was to develop and demonstrate technologies which improve 
the etfoctiveness and efficiency of water quality monitoring programs. 


rarr 


LANDSAT SATELUTE CHARACTERISTIC S 


The ^andsat aystem was developed by the National Aeronautics and 
Space Administration (NASA) to help meet the increasing demand for man 
to manage the earth's limited natural resources. Data gathered by the 
Landsat has been applied to studies in the fields of geology, cartography » 
geography, land management, forestry, hydrology, and many others. 

This series of satellites began with the launch of Earth Resources 
Technology Satellite I (ERTS) In July 1972. This experimental satellite 
proved the applicability of monitoring the earth'.*; surface from space, 
and led to the launch of ERTS 2 in January 1975. It was after the 
successful launch of ERTS 2 that the new name Land Satellite (Landsat) 
was adopted. The new name distinguishes these satellites from the Seasat 
series of earth observation satellites. 

i.andsat 1 was turned off in January 1978. landsat 2 continues to 
operate and landsat 3 was launched in March of 1978. Plans are currently 
being developed for the fourth landsat with a scheduled 1981 launch. 
Imagery used in this study was obtained by landsat 1 or landsat 2. 

landsat satellites arc launched into sun>synchronous near-polar 
orbits at an altitude of approximately 900 km (540 mi) (table 1). 

This type of orbit ensures repeatable sun-illumination conditions for 
any particular date from year to year. 

The satellites cross the equator every 103 minutes thus completing 
14 orbits in 24 hours. Therefore, the next westward track of data for 
any orbit is acquired at the same sun time the following day (figure 1). 
The earth rotates 2,760 km (1,650 mi) under the satellite at the equator 
during each orbit. The coverage width of each orbit pass is 185 km 
(115 mi) and the distance between adjacent orbits at the equator is 
159 km (95 mi). Complete earth coverage is, therefore, completed by 
each satellite every 18 days, landsat 2 was launched so that its orbit 
follows landsat 1 by 9 days, landsat 2 and 3 also provide 9-day coverage. 

The instrumentation of landsats 1 and 2 consists of two imaging 
systems, the multispectral scanner (MSS) and the return beam vidicon 
(RBV), Also on board are the data collection system (DCS) receiver and 
transmitter, and two wide band video tape recorders (figure 2). Only 
data from the MSS system was used in this study and need be considered. 

The MSS is a line-scanning radiometer which collects data by creating 
images of the earth's surface in four spectral bands simultaneously. 
Radiation coming from the surface of the earth and its atmosphere is 
recorded as an analog signal which is converted to values of from 0-63. 

The numbers represent brightness values (BV), the amount of electromagnetii 
energy reflected from an area on the earth's surface in one wave length 
band. 


The MSS scans the earth's surface from west to east (figure 3). 
Twenty-four detectors are used to record six lines of data (figure 4) in 
each of the four wave length bands (table 2), 
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VIOICON CAMERAS (3) 


Figui'e 2; Landsat configuration. From NASA landsat Data Uiera Handbook 
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During •* ulgnal vi$ sampled every OtOS microseconds. For 

iMch hiiiid, .ippi’oxiimiiely 1, 300 samples are Lakeii along a IB5 km line 
(ligure 4). Tints, the inslant.incous field ol view (IFOV) of 79 m by 
79 III moves about 56 m on the ground between each sample. The individual 
radiation measnremeiits mitst he arranged on an image in a manner that 
preserves spatial relationships. Thiis, the ineasnrements arc assigned 
dimensions ol .56 m by 79 m so that geometric distorlLons are not 
introduced. The 56 m by 79 m area is called a Landsat picture element 
or pixel. )\ 

Landsat MSS imagery is placed in the public domain and is available 
as either pliotographic products or computer compatible magnetic tape. 

For the ii.ser to locate the are.i of bis interest, the continuous 
image of the MSS has been divided along the orbit path (north to south) 
into sections eipial to the east-west width of the MSS scan, 185 km 
(115 ml). This division is always made as near to the same location 
as possible, thus creating nominal scenes of Landsat data. These scenes 
are assigned a unique identifying number corresponding to the orbit path 
and the east-to-west row of scenes (figure 5). 

Pliotographic products available include black and white prints of 
individual Landsat hands. These products cover one nominal scene of 
imagery 185 km by 185 km and are available in a variety of scales. Also 
available are false color infrared composites of selected scenes. These 
products utilise band 4, 5, and 7 to create photographically the false 
color image. 

Landsat computer compatible tapes (CCT) are available in one-tape, 
1,600 bit-per-inch (bpi) or in two-tape 800 bpi format. 


Nowinul scenes foi' Londsat ijnagei'y 





tkn nesske valley keoioh study area 


This study focuses on tho gcester Tennessee Vslley tegion of the 
southeastern United States, The drainage basin of the Tennessee River 
<Micom|)(iSKes n land-locked area of nhont 1Q6»000 sq. km (40,910 sq. mi) 
ill the soutlieiiuiern United Status, inuLuding |tartH of Tennussue, Alahnma, 
Horth CaroUna, Virginia, Georgia, Kentucky, and Hississippl. 

The Tennessee River system has almost 61,200 km (38,OOC mi) of 
streams and rivers, bakes and reservoirs have over 2,600 sq, km 
(1,015 sq, mi) of surface water and more than 17,700 km (11,000 ml) 
of shoreline. The w,atershcd is characterissed by rugged mountains and 
green forests in the eastern portion of the Valley and rolling hills, 
open fields, and woodlands in the west. From Mount Mitchell, North 
Carolina, in the east, to Paducah, Kentucky, in the west, the topography 
ranges from 2,037 m (6,684 ft) to 90 m (300 ft) above sea level, 

The Tennessee River Valley, one of the wettest regions of the 
United States, receives about 132 cm (52 in) of rainfall a year, March 
is usually the wettest month, and September or October, the driest. 

The climate is mild and humid with an average mean air temperature of 
about 15®C (59®F) and an average mouthly humidity of 66 to 84 percent. 

About $9 percent of the Tennessee River basin is forested; about 
38 percent is open land and pasture; and 2 percent is covered % water. 
Approximately 4 million people live in the Tennessee River watershed. 

Over 85 percent of this population resides in towns and cities or in 
the nonfarin rural areas surrounding population centers; fewer than 15 
purcunt now live on farms. About 40 percent of the population lives in 
six population vtMiiers! Asheville, NorLh Carolina; Hi'istol-Johnson City- 
Kingston, Tennessee; Knoxville-Oak Ridge, Tennessee; Chattanooga, 
Tennessee; Huntsvllle-Decatur, Alabama; and Florence-Sheffield-Tuscumbia, 
Alabama . 

The reservoirs examined in this study are located in parts of 
Tennessee, Alabama, Georgia, Kentucky, and North Carolina (figure 6), 
Superimposed on this figure of the greater Tennessee Valley region are 
the approximate ground area coverages of t,a»dsat scenes corresponding 
to Path 20, Row 35 8 36 (figure 5) from which information relating to 
reservoir trophic status was extracted. This region Was selected because 
of the many ivservoirs and the availability of ground truth water quality 
data and Landsat imagery. 

Ground truth data for the reservoirs incorporated in this study were 
collected by the U.S. Environmental Protection Agency (EPA) during their 
National Eutrophication Survey (NES) of 1973. Tlic study reservoirs are 
liited In table 3 with information of the morphometry and hydrology of 
selected reservoirs, NES-sampled reservoira tra the first 35 iitted; 
the last 14 listed reservoirs were outside the scope of the NES, but 
are included becauie they are of Intereat In using Landsat to sssest 
their trophic state, 







lY AND HYDROLOGY OF 35 RESERVOIRS IN THE GREATER TENNESSEE VALLEY REGION 
SAMPLED DURING THE NATIONAL EUTROPHICATION SURVEY OF 1973 
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MORPHOMETRY AND HYDROLOGY OF 35 RESERVOIRS IN THE GREATER TENNESSEE VALLEY REGION 
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GROUN D TKUTH WATER QUALITY DATA COLLECTION 


Ali lh«» 35 NKS fc'servoirs were sampled three times (spring, SMiwiyer, 
and fall) during the 1973 ccl«ndar year by pontoon*equipped helicopter- 
borne sampling teams. The helicopters were equipped with in situ sensors 
for the measurement of conductivity, temperature, optical transmissivity, 
dissolved oxygen, pit, and water depth. Samples for algal identification, 
chlorophyll a, and nutrients were collected using a submersible pump. 
Additional equipment included an echo sounder, 30 cm Secchi disc, and 
water sampling equipment, Specific details of sample collection and 
methods of laboratory analyses are found in National Eutrophication 
Survey Methods (1975). 

Most of the reservoirs surveyed were chosen on the basis of actual 
or potential eutrophication problems, with the result that this inves- 
tigation does not universally represent the normal distribution of 
reservoirs in the Tennessee Valley, with respect to trophic state, and 
is biased or weighted toward reservoirs that are often referred to as 
"ciitrophic." In general only reservoirs of one km^ or larger and mean 
hydraulic retention times of at least two weeks were considered; however, 
these selection criteria were waived for reservoirs of pattieiilar 
intere.st. Sampling .sites in each reservoir were selected primarily to 
attempt to deline the char.scter of the reservoir as a whole rather than 
specific areas or embayments of the reservoirs and chosen to reflect 
the deepest portion of each major basin in a reservoir. The number of 
sampling sites varied for different reservoirs, ranging from one 
(Jimaluska Reservoir) to seventeen (Kentucky Reservoir). 

Data mve made available through the Environmental Protection 
Agency's water quality data storage and retrieval system (STORET). 


EUTROPHICATION AND TROPHIC STATE INDICATORS 


Vollouweidcr (1968) has defined eutrophication of water bodies as 
"• . .their enrichment in nntrient.<; and the ensuing progressive deteriO'^ 
ration of their ()uality, especially lakes, due to the luxuriant; growth 
of plants with its repercussions on the overall metabolism of the waters 
affected. . 

Most lakes and reservoirs originate as water bodies possessing 
relatively low concentrations of nutrients and generally low levels of 
productivity. As the water body ages, Inflows carry sedimont which 
decreases a water body's depth, and nutrients, which stiiimlate produc- 
tivity and further increase the sedimentation rate. Floral and faunal 
changes occur. Algal blooms become more coiwnon, rooted aquatic species 
increase, and desirable game fish may be replaced by rough fish. 

Naturally, this eutrophication process is very slow with Llio norm.il 
life span of a lake being on the order ot several hundred yo.Trs, How- 
ever, man'.s pr.'ictices relating to the disposition of nutniclp«ti .sewage, 
industrial wastes, and land use cause large nutrient loads on many water 
bodies. This results in a rapid aging of the lakes and reservoirs and 
makes the water bodies less attractive to users and, more importantly, 
shortening the life of the lake or reservoir. 

Many different physical, biological, and chemical cliaracteristic.s 
are needed to adequately describe a water body's trophic state. There 
are many indic.itors of trophic state, each with its merits and short- 
comings and many opinions regarding which indicator(s) should he used 
in classifying lakes and reservoirs. Table 6 lists some of these common 
indicators or Indices. These many indicators of trophic state reflect 
the multidimensional problem of classifying water bodies and defining 
its trophic state condition. 

The historical aspects and sematics associ.ated with the words 
"eutrophication," "oligotrophic," "mesotrophic," and "eutrophic" are 
found in Weber (1907), Naumann (1919, 1931), Thienemann (1918), Rodtie 
(1969), Hutchinson (1967, 1973), Becton and Edmondson (1972) and Edmonds 
(1974). 
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WATER QUALITY DATA ANALYSIS 


Ffoir the very beginning ot this study it vns apparent that there was 
a need to provide a more realistic assignment of a water body's trophic 
condition than the nebulous and overlapping categorissations of oligo- 
trophiCi mesotrophic, cutrophic, or hypcreutrophic. In a large group of 
lakes I trophic condition is a continuum with no sharp classes as suggested 
by these classical groupings. 

A relative » numerical tropnic state index, based on ground truth 
survey data was needed to better define trophic condition (i.e*) 
quantitative rather than qualitative). 

Description of a lake or reservoir's trophic state requires consid* 
eratioii of several difficult physical, biological, and chemical 
characteristics) and as such trophic state cannot be directly measured 
in the field. 

Because of the very nature of the multidimensional concept of 
trophic state, trophic classification lends itself well to two multi'* 
variate statistical techniques, cluster analysis and principal component 
analysis. Brezobik and Shannon (1971) were among the first to apply 
multivariate techniques in their classification of 55 lakes in Florida, 

Thirty-five resjhrvoirs, sampled by the NF.S in 1973 in the greater 
Tennessee Valley reyjion, were selected for water quality data analysis. 

A careful examination of the water quality characteristics measured by 
the NES and a review of pertinent literature (Boland, 1976; Brezonik and 
Shannon, 1971; Shannon and Brezonik 1972a, 1972b; Carlson, 1977; Hooper, 
1969; lueschow, 1970; EPA, 1974; and Rast and bee, 1978) resulted in the 
selection of six indicators of trophic states conductivity, umhos/cm; 
chlorophyll a, pg/1; total phosphorous, mg/l, total ocganic nitrogen, 
mg/l; algal assay yield, dry-weight in mg/l', and Secchi disc transparency, 
inches. So that all indicators would contribute to the trophic state in 
a positive sense (i.e,, increasing value of indicator being associated 
with increasing eutrophication) inverse values for Secchi disc depth were 
used in the data analyses and in development of the trophic state iiulices 
for each reservoir. Annual mean values were used in the aiialyses. A 
lack of normality in the data necessitated a natural log transfonnation 
of the mean values prior to the data analyses , The annual mean values 
for the six trophic state indicators for each of the 35 reservoirs are 
given in table 5 with descriptive statistics. Table 6 is a correlation 
matrix of these six trophic indicators iu which the coefficients were 
determined using natural log transformed data for the 35 NES sampled 
reservoirs. 


Cluslcr An n ly.sl.s ■ 

Cluster anaVysis was used to find groupings of reservoirs with 
similar trophic states, based on the six indicators given in ial)le 5. 
This study eniployed the NT-SYS program for complete linkage cinstering 
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(Rohlfi Klffhpaugh, And Ktrk, 1971) iiiing the Euclidiiii dittanct to txMint 
the 35 reiorvoira aamplad during the 1973 NE8 for natural groupinga. The 
AM'lliud Ik hIko hitown na the Ciirtheat neighbor method. An excellent revitm 
of the method ia given in Sneath and Sokai (1973). 

The reiulta of the cluatec analyaia are ahown aa a dendogram in 
figure 7. The abaciaaa ia the Euclidian diatance. On the ordinate axia» 
the auihora have attempted to rank the reaervoira according to trophic 
atate by uaing atom rotation. In general, the trophic atatiia increaaea 
along the ordinate in a upward direction. 

What we have atteiq>ted to illuatrate with figure 7 ia that, depending 
upon the criteria uaed to define "a cluater," theru are certainly more 
than the three claaaic atatea (eutrophic, meaotrophic, and oligotropbic) 
of reaervoir eutrophication and that aa the number of reaervoira being 
conaidered increaaea a trophic "continuum" developa. There ia aome diffi- 
culty in reconciling the eight cluatera (A, B, C, li, E, F, G, H) with the 
three claaaic trophic atatea. Uaing KES aaaersmenta aa a guide: Cluater 
A may be considered as a very hype rent rophic lake; clusters B, and D 
may be considered as a mixture of hype reut rophic and eutrophic reservoirs; 
clusters E, F, and G may be considered as a mixture of eutrophic and meso- 
trophic reservoirs; and cluster H consists of reservoirs characterised as 
both mesotrcphic and oligotropbic. 


Principal Component Analysis 

Principal component analysis is an ordination technique which in 
effect reduces the concept of trophic state from six indicators to a 
single index. In general terms, principal component analysia reduces 
the dimensions of a concept (eutrophication) by expressing the original 
observations (six trophic state indicators) in fewer terma (trophic state 
index). The first principal component (PCI) of a set of variables is the 
linear combination of the variables which explains the maximum variance 
in the original data. Ho absolute physical meaning can be placed on the 
trophic state index (PCI) Values, but the values are felt to have relative 
meaning and are a quantitative expression of the trophic condition of the 
35 reservoirs sampled by the 1973 NES in the greater Tenneasee Valley 
region. The NT-SYS system was used to perform the principal component 
analysia, (Rohlf, Kishpaugh, and Kirk, 1971). Detailed descriptions of 
the theoretical and computational aspects of principal coaq>onenta are 
found in Hotelling (1933a, 1933b, 1936), Anderson (1958), and Horrison 
(1967). 

The normalised eigenvalues and eigenvector numbers are given in 
table 7. The first component (eigenvector 1) accounts for about 67 
percent of the variation in the data. Correlation coefficients between 
the principal components and the six natural leg transformed trophic 
indicators (table A) show that the first principal component (PCI) ia 
highly correlated with each of the trophic indicators. From the results 
shewn in tables 7 and 8 it was concluded that the first principal 
component ia indeed indicative of each reservoir's relative position 
on a multivariate trophic scale. 
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Tttble 7 


A 

NormoXiKQd Eigenvalues and Elgenvectore 


Eigenvector 

Number 

Eigenvalue 

Variance 

Cummulativc 
Variance (X) 

1 

4.012 

66.87 

66.87 

2 

0.921 

15.35 

82.22 

3 

0.493 

8.21 

90.43 

4 

0.270 

4.50 

94.93 

5 

0.218 

3.63 

96.56 

6 

0.086 

1.44 

100,00 


^I'he principal component analysis was performed using naturol log 
transformations of mean data for six trophic state indicators 
(Inverse Secchlt chlorophyll a, algal assay » conductivity, 
total organic nitrogen, and total phosphorus) for the 35 
reservoirs sampled during the 1973 NES in the greater 
'I’enneasee Valley region. 


•i’abiQ 0 

Covtvltttton Co«rriv‘l«*nC8 of 'fpophic 
State Indicators and t'rliKipal Comnonenls 


Principal Coinponont 


Inverse Socchl 

0.879 

0.2U 

“0.161 

-0.268 

0,312 

-0. 

Conductivity 

0.755 

0.366 

0.688 

0,231 

0,051 

-0. 

Total Organic 
Mtroaan 

0.758 

-0.672 

0.313 

0.277 

-0. 166 

-0, 

Total Phosphorus 

0.957 

0.07fi 

-0.U6 

0,033 

-0.037 

0. 

Chlorophyll ^ 

0.682 

-0.66S 

-0.161 

0.275 

0.117 

-0, 

Algal Assay 

0.8/.5 

0.310 

-0.312 

0.061 

-0.276 

-0. 


1973 in fell® greater TunnesRao Valley nop lnv\, 
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Tae reservoir with the lowest PCI value, Burton, is rated ss having 
the lowest trophic state of those studied. Trophic state increases in 
the positive direction with Reelfoot exhibiting the highest trophic ! 
state of the 35 water bodies studied. As further evidence that the PCI 
values have real meaning, normalized mean composite rank (NHR) indices 
were calculated and compared as shown in table 9. The HHR's were cal'’ 
culated by ranking each of the six indicators for each reservoir frPR 1 
through 35 and then calculating the average rank for the six indicators 
at each reservoir. Finally, these averages were normalized. It is quite 
evident that the NMR and PCI Values results in essentially identical 
rankings for the 35 reservoirs. This ranking approach has been used in 
other similar studies (EPA, 1974; Lueschow, 1970; and Piwoni and Lee, 

1975) . 

In conclusion, the PCI's presented in table 9 represent an assessment 
of each reservoir's relative trophic state and were used as trophic state 
indices to evaluate Landsat-reservoir eutrophication relationships. 


•rtible 9 


Principal Component Value and Norinallxed Moan Rank 
Index lor $5 Koncrvolrs Sampled Curing t!io National Eutrophication 
Survey of 1973 In the Greater Tennessee Valley Ueglon 


Reservoir 

Name 

Identification 

Number 

PC 1 
Value 

NMll 

Index 

Positlo 

PCl/NMR 

Burton 

GA 6 

-1.441 

-1.720 

1/1 

Blue Ridge 

GA 4 

-1.274 

-1.679 

2/2 

Snnteetlah 

NC 27 

-1.216 

-1.516 

3/3 

Chatuge 

GA 7 

-1.051 

-1.354 

4/4 

Dele Uollou 

KY 12 

-0.970 

-1.130 

5/7 

Fontana 

NC 14 

-0.960 

-1.231 

6/6 

Nottely 

GA 22 

-0.957 

-1.242 

7/5 

Hlwassee 

NC 18 

-0.767 

-0.967 

8/9 

Cumberland 

KY U 

-0.612 

-1.058 

9/8 

Sidney Lanier 

GA 28 

-0.524 

-0.702 

10/10 

South Hols ton 

TN 29 

-0.503 

-0.590 

il/U 

Tims Ford 

TN 30 

-0.368 

-0.570 

12/12 

Alin toons 

GA 1 

-0.266 

-0.041 

13/16 

Woods 

TN 35 

’“0*46 A 

-0.315 

14/14 

Great Falls 

TN 16 

-0.264 

-0.132 

15/15 

Chickamouga 

TN 10 

-0.153 

-0.397 

16/13 

Watts Bar 

TN 32 

-0.036 

-0.163 

■' 17/21 

Douglas 

TN 13 

0.075 

0.051 

18/17 

Wilson 

AL 34 

0.187 

0.092 

19/18 

Pickwick 

AL 25 

0.191 

0.122 

20/19 

Junaluska 

NC 19 

0.258 

0.153 

21/20 

Nickajack 

TN 21 

0,291 

0.214 

22/22 

Barren 

KY 3 

0.390 

0.743 

23/25 

Boone 

TN 5 

0.441 

0.712 

24/24 

GuntersviUe 

AL 17 

0.478 

0.641 

25/23 

Fort Loudoun 

TN 15 

0.530 

0.855 

26/27 

Old Hickory 

TN 23 

0.581 

0.753 

27/26 

Kentucky 

KY 20 

0.612 

0.967 

28/21 

Percy Priest 

TN 24 

0.694 

1.109 

29/2<i 

Cherokee 

TN 9 

0.714 

1.160 

30/3( 

Cheatham 

TN 8 

0,908 

1.333 

31/34 

Barkley 

TN 2 

0.929 

1.252 

32/33 

Weis 

AL 33 

1.004 

1.201 

33/31 

Waterville 

NC 31 

1.067 

1.252 

34/3 ‘ 

Reelfoot 

TN 26 

2.278 

1.873 

35/3' 


li 




italic 


-25 


Cion 
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*2 

3 

n 

'6 

/5 

•9 

/8 

/lO 

/ll 

/12 

/16 

/14 

715 

713 

/2l 

m 
/1 8 
719 
./20 
722 

725 
724 
723 
727 

726 
1/28 
729 
)/30 
1/34 

732 
1/31 

733 
)/35 


L ANDSAT DATA EXTRACTION 


Imagery Selection and Manipulntioil 

The extraction of Landsat data |or each reservoir required the 
selection of appropriate dates of imagery, defining the reservoir 
location on each scene, and then computing the necessary statistics 
of the raw landsat data for each reservoir. The computerized data 
processing portion of this work was performed at the University of 
Wisconsin, by the Environmental Monitoring and Data Acquisition Group. 

Landsat scenes were selected to meet the following criteria: 
maximum number of reservoirs on the minimum number of scenes; minimum 
cloud cover; date close to time of sampling; and good quality imagery 
(figures 8-12). These criteria were difficult to meet and led to poor 
correlation between dates of Landsat coverage and NES water sampling 
as no acceptable scenes were found for 1973. On the majority of dates, 
cloud cover was excessive (greater than 10 percent). Even using scenes 
with up to 50 percent cloud cover did not provide adequate coverage in 
1973. This problem is shared with others, but the Smokey Mountains are 
in fact deserving of their name. 

With the launching of Landsat 2 and the resultant 9-day coverage 
pattern, adequat<> Landsat coverage has been recciv(si each year since 
1975. 


The search for Landsat imagery was performed at the EROS Browse 
File, operated by the TVA Mapping Services Branch, from l6mm microfilm. 
Computer compatible magnetic tape and black and white photographic 
products were obtained for the selected scenes. Four sets of tapes 
were obtained through the courtesy of Oak Ridge National Laboratory and 
one was obtained directly from the EROS Data Center in Sioux Falls, 
South Dakota. The scenes selected and used in this study are listed in 
table 10. The reservoirs which were examined in this study are listed 
In tables 11 and 12. 

The first step in locating each reservoir on a Landsat scene was 
to construct a grid overlay which could be used to read the Landsat 
coordinates of any point on the scene. This grid overlay or **pixel 
counter" was drawn on mylar to fit each frame of Landsat data. It was 
then possible to find each reservoir on the photographic image and 
obtain the range of scan lines and elements that would cover that 
reservoir. This method of determining Landsat coordinates was accurate 
to approximately 50 pixels. 

These Landsat coordinates were used to generate a line printer 
character plot for each reservoir (figure 13). These plots were for 
band 7 caw reflectance values. Values of 0-7 were left blank and were 
used to refine the polygon describing each reservoir, eliminate shadow 
effects, and to verify the reflectance level for that reservoir. 
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figuie 8 Reproduction of LANDSAT Multispectral Scanner data, 

Band 7 (0 8 1.1 micrometers) Frame 1084 15431 (15 October 1972) 
Identified are classified lakes. Lakes underlined were sampled as part of 
National Futrophication Study m 1973 
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Figure 9 Reproduction of LANDSAT Multispectral Scunner data. 

Band 7 (0 8 11 micrometers) Frame 1822-15315 (23 October 1974) 
Identified are classified lakes Lakes underlined were sampled as part of 
National Lutrophication Study in 1973 
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Figure 10 Reproduction ot LANDSAT Multispectral Scanner data. 

Band 7 (0 8 1 1 micrometers) Frame 1822 15322 (23 October 1974) 
Identified are classified lakes. Lakes underlined were sampled as part of 
National Eutrophication Study in 1973 
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Figure 11 Reproduction of LANDSAT Multispeclral Scanner data. 

Band 7 (0 811 micrometers) Frame 1948 15264 (26 February 1975) 
Identified are classified lakes. Lakes underlined were sampled as part of 
National Eutrophication Study in 1973. 
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TabU 10 

LANDSAT MSS FRAMES 


Frasie 

Nujsber 

Date 

Area 

Nusiber of 

Lakes 

1084-15431 

15 October 1972 

Northeastern Tennessee 
Western North Carolina 

12*» 

1822-15315 

23 October 1974 

Northeastern Tennessee 
Western North Carolina 

14" 

1822-15322 

23 October 1974 

Southeastern Tennessee 
Western North Carolina 
North Georgia 

IS 

1048-15264 

26 February 1975 

Northeastern Tennessee 
Western North Carolina 

is" 

2224-15303^ 

3 September 1975 

Northeastern Tennessee 
Western North Ciirolins 

is" 


Jlandsat 2 

‘’lO reservoirs are convnon to all four scenes. EPA data was collected on 
6 of these reservoirs. 
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Table 11 


DATES OF tANDSAT DATA FOR SAMPLED RESERVOIRS 


Reservoir Name* 
and Stale 


Allntoona 

GA 

1 

Rlue Ridge 

GA 

4 

Boone 

TN 

5 

Burton 

GA 

6 

Chatuge 

GA 

7 

Cherokee 

TN 

9 

Chlckamauga 

TN 

lO 

Douglas 

TN 

13 

Fontana 

NC 

lA 

Fort Loudoun 

TN 

15 

Hiwasspe 

NC 

18 

Junaluska 

NC 

19 

Nottely 

GA 

22 

Santeetlah 

NC 

27 

Sidney Lanier 

GA 

28 

Watcrville 

NC 

31 

Watts Bar 

TN 

32 




in 

m 






Cv 







a 

» 

1 

1 

iJ 

4J 




s 

V 

t 

w 

wn 
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iN 



X 

X 

X ■ X 

X 
X 

X X X X 

X 

X X X X 

X X X X 

X X X X 

X 

X X 
X 

X X X X 

X 

X 

X X X X 


*Reservoirs sampled by the Environmental 
Protection Agency in 1973. 
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Tabic 12 

nATES OF bAKDSAT DATA FOR NONSAMPLED RESERVOIRS 


Reservoir Name*" 
and State 




in 

in 





«7v 

o\ 

a\ 



r-d 



!< 

• 

» 

* 

U 

U 


a 

u 

u 


Oi 

o 

p 


<8 

in 

m 




esi 
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Carters 

GA 

1 

Chilliowee 

TN 

2 

Hartwell 

GA 

3 

laurel 

KY 

4 

Helton Hill 

TN 

5 

Nantahala 

NC 

6 

Nolichucky 

TN 

7 

Norris 

TN 

8 

Parksviile 

TN 

9 

Rabun 

GA 

10 

Thorpe 

NC 

11 

Toxsway 

NC 

12 

Tuckasegee 

NC 

13 

Tugaloo 

GA 

14 


••'‘Reservoirs included in this study which were 
not siimpled by the Environraental Protection 
Agency in 1973, 
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Figure 13; Computer printout 
of Landsnt band 7 density 
slice. Shown is s portion 
of Cherokee Reservoir from 
scene 1822-15315, The area 
labeled 01 is a NES sample 
Bite. 
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To determine which picture elements were actually in the water a 
band 7 deniity slice was made. Only those picture elements with band 7 
values of 0>3 were used as representative of wai,er. This relatively 
low value as compared with those used in other studies (Scarpace 1$76» 
Boland 1976) was necessary for several reasons. This region of the 
southeast is very rugged with ridges and valleys, which cause signifi- 
cant shadows. These prominent shadows exhibit low reflectance readings 
that may be confused with water. This effect is reduced using a 0-3 
density slice. Another advantage to using the 0-3 band 7 density slice 
was the elimination of shallow water bottom and shoreline effects. This 
also resulted in improved correlation between the ground truth water 
quality data and the MSS data because all NES water samples were col- 
lected in deeper water. From inspection of the character plots, all 
pare water pixels showed reflectance values of 3 or less in hand 7. 


Generation of Statistics 


All picture elements falling within the polygon which generally 
outlined each lake were extracted from the raw data tape, Of these 
picture elements, only those with band 7 values between 0-3 were con- 
sidered in the statistical analysis. Fourteen statistics were deter- 
mined for each reservoir. These are the mean raw reflectance values 
for each band (MBA, MBS, HB6, MB7), the variance in each band (VBA, VB5, 
VB6, VB7), the ratio between band A and 5 (RBA), band 5 and 6 (RB5) and 
band 6 and 7 (RB6), and the variance of these three ratios (VRBA, VRB5, 
VRB6). The ratios were computed using the following formula: 

RBi a ^ ARCTAN (Bi/Bi+1) 

where 

RBi a ratio between band i and i-fl 
Bi a raw reflectance of band i 
Bi+1 a raw reflectance of band i+1 

This pic'duces a range in the value of each ratio from 0 to 255 where a 
change of one unit is always equivalent to the same change in direction 
in spectral space. 

Estimates of three trophic state indicators (Secchi disc depth, 
conductivity, and total phosphorus) and trophic state indices (PCI) 
using Landr It and NES ground truth data are demonstrated in the remainder 
of this report. 

In an effort to reduce the amount of data presented in the text, 
attention has been focused pn frames 1822-15315 and 15322. Prior to 
processing, these two scenes were in actuality one continuous frame 
recorded on the same orbit. For this reason scenes 1822-15315 and 
15322 are treated as one. Regression models and data for frames 
222A- 15303, 19A8-1526A and 108A-15A31 are presented in Appendix A. 


Bi 


RELATIONSHIPS BETWEEN TROPHIC STATE INDICATOR AND lANPSAT IMAGERY 


Lartdsat cannot, ditfcctly inoasui'o cheinical indicntor^ of water quality 
but its areal and spectral resolution porinit the detection of pUenoiticna 
indirectly related to eutrophication. In each frame data were extracted 
for each of the four bands for the NES-aampIcd reservoirs . The mean and 
the variance of the reflectance values for each band, and the ratio and 
the variance of the ratio between bands were calculated. A correlation 
analysis was performed using this Iiaiuisat data and the 15 NES reservoirs 
in frames 1822-15315 and 15322. Correlations between the Landsat data 
anu trophic state index (PCI), Secchi di.sc deptli, conduetivity, and total 
phosphorus are found in table 13. Several correlation coefficient pairs 
(e.gi, band 6 and Secchi disc depth) exhibit high correlations inferring 
that relationships do exist between L.’indsat reflectance values and water 
quality characteristics. 

Data analysis and regression models were developed using the Statis- 
tical Analysis System (SAS) (Barr, Goodnight, Sail, and Helwlg, 1976); 
the maxinuun R“ Improvement technique oi the stepwise multiple regressloii 
procedure was used to develop multilinear regression models. As such 
numerous models were developed in this investigation. Criteria used in 
the selection of the "best" models included (the magnitude of the 
square of tlie multiple correlation coefficient) and the F-statistic. All 
regression coefficients were required to be significant at the 0.05 level. 


Secchi Disc Depth Estimation 

The best regression model, as measured by the square of the multiple 
correlation coefficient (R^) and the F-statistic, for estimating Secchi 
disc depth is: 

log^ (Secchi) = 6.282 + 0.162 MBS - 0.598 MB6 + 0.006 VRB4 

The model accounts for about 96 percent of the variance about the mean, 
table U. The observed and predicted Secchi disc depth values for the 
15 reservoirs ace given in table 15. 

Althougli the models presented hero are purported to estimiite Secelii 
disc depth, couductivif$, and total pho.sphorusV it must be remembered that 
the ground truth data collected and Landsat overflight are not concurrent, 
being Separated in time by about one year. This "nonconcurrence" limits 
the reliability of the iaodeis and precludes more precise estimates. Cau- 
tion must be exercised In .assuming that the models are applicable to other 
reservoirs or even to the same 15 reservoirs on a different date, Reservoirs 
are by their very nature dynamic. In a period of a few days or weeks, 
their appearaiijce c.an change significantly due to algal blooms, turbidity 
plumes caused '^y heavy rains, and man induced changes in reservoir volume. 

In addition to/reservoir dynamics are variations caused by atmospheric 
conditions anji solar angle. As such the models .‘should be tceated a.s 
"ball park" rather than accurate estimators. In fact it is remarkable 
that estimates are as good as they are and further emphasises the fact 
that relationships do indeed exist between water quality cUaracteristics 
and Landsat Imagery. 
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■Jdble 13 

>i 

CoiToXntJons Between Ground Truth Wnter Quality Bata 
and 1.AN03AT Dnta Cor 15 Reservoirs in Vrames 
18?.2«15315 and 15322 


Band A (MBA) 

• 

PC 1 
O.AAi 

Band 5 (MB5) 


0.618 

Band 6 (HB6) 


0.808 

Band 7 (MB7) 


0.72A 

Band A/Bnnd 5 (RBA) 

1 

-0.818 

Band 5/Band 6 (RB5) 

1 

1 

1 

-0.066 

Band, 6/Band 7 (RB6) 

0.386 

Visrianco Band A (VBA) 


0.130 

Vaviancc Band 5 (VB5) 


0.128 

Variance Band 6 (VB6) 


0.118 

Variance Band V (VB7) 


-0.77A 

Variance Band A /Band 5 

(VKBA) 

-0.573 

Variance Band 5/Band 6 (VRB5) 

-0.686 

Variance Band 6/Band 7 

(VRB6) 

-0,730 


lOBft 

Secchi 

Gonductivitv 

Phosphorus 

-0.5A7 

0.557 

0.A83 

-0.698 

0,717 

0.6A3 

-0.918 

0,863 

0.783 

-0.853 

0.697 

0.6A6 

0.831 

-0.858 

-0.821 

0.1 2A 

0.028 

0.0A3 

-0.313 

0.A08 

0.A06 

-0,366 

0,287 

0.131 

-0.36A 

0.283 

0.128 

“0.355 

0.273 

0.118 

0 . 79 , A 

-0.8A9 

-0.732 

0.582 

-0.595 

-0.53A 

0.705 

-0.722 

-0.727 

0.771 

-0.763 

-0.768 
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Analysia of Vari.iticc 

Source DF Sum ol Squai*es"*~**^^* Mean Square 

Regression 3 3.481 1.160 103.98 


Residual 

Total 

II 

14 


0.123 

3.604 

0.011 

(R2 a 0.966) 

®This analysis 


nerformed using Lands.nt 

data from frames 

1822-15315 

and ir>322. 
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Table 15 





.SECCIII DISC DEPTH RESIDUALS® 





Observed 

Predi eted 





Seech i Disc 

Secchi Disc 

Residual 

Reservoir 

IdentiFication 

1 Deiuh 

DeptI* Observed-Predicted 

Name 


Number 

(inches) 

(inches) 

(inches) 

Allatoonn 


OA 1 

56.8 

53.1 

3.7 

111 no Ridge 


GA 4 

105.1 

109.2 

-4.1 

llurLon 


GA 6 

136.1 

127.0 

9.1 

Cliatuge 


GA 7 

117.2 

109.2 

8.0 

Cherokee 


TN 9 

51,3 

48.1 

3 . 2 

Chtckaniauga 


TN 10 

36.2 

35.8 

0.4 

Douglas 


TN 13 

57.2 

62.8 

-5.6 

Fontana 


NC 14 

107,4 

104.6 

2.8 

Fort Loudoun 


TN 15 

34.4 

31.8 

2.6 

lliwassee 


NG 18 

79.4 

90.5 

-11. 1 

Junaiuska 


NC 19 

38.0 

39 . 1 

-1 . 1 

NoLLoly 


GA 22 

94.3 

96.7 

-2.4 

Santeetlah 


NC 2? 

133.0 

11 1.0 

22.6 

Sidney Lanier 


GA 28 

103.6 

120,0 

-16,4 

W.1LLS B.U* 


TN 32 

37.8 

42,2 

-4,4 

“'This ana lysis 

was 

performed 

using Landsat 

data from frames 

1822M5315 and 
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Conductivity Estimation 

The best regression model found for prediction of conductivity is; 

log^ (conductivity) = 6.418 + 1.580 MB7 - 6.848 VB7 

The model accounts for about 85 percent of the variance about the meant 
table 16. The observed and predicted conductivity values for the 15i 
reservoirs are given in table 17. 


Table 16 

ANALYSIS OF VARIANCE - CONDUCTIVITY® 




Analysis of Variance 

Source 

DF 

Sum of Squares 

Mean Square 

F 

Regression 

2 

13.115 

6.557 

33.14 

Residual 

12 

2.374 

0.198 


Total 

14 

15.489 


(R2 « 0.847) 

®This analysis 

Was performed using Landsat 

data from frames 

1822-15315 

and \5322. 
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TABLE 17 
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CONDUCTIVITY RESIDUALS® 




Observed 

Predicted 

Residual 

Reservoir 

Identification Conductivity Conductivity Observed-Predicted 

Name 

Number 

(pmhos/cm) 

(pmhos/cm) 

(pmhos/cm) 

Allatoona 

GA 1 

40 

65 

-25 

Blue Ridge 

GA 4 

16 

31 

-15 

Burton 

GA 6 

17 

21 

-4 

Chatuge 

GA 7 

21 

33 

-12 

Cherokee 

TN 9 

260 

124 

136 

Chickamauga 

TN 10 

158 

213 

-55 

Douglas 

TN 13 

184 

93 

91 

Fontana 

NC 14 

27 

27 

0 

Fort Loudoun 

TN 15 

210 

203 

7 

Hivassee 

NC 18 

32 

24 

8 

Junaluska 

NC 19 

106 

137 

•31 

Nottely 

GA 22 

25 

17 

8 

Saatcetlah 

NC 27 

17 

13 

4 

Sidney Lanier 

GA 28 

35 

35 

0 

Watts Bar 

TN 32 

176 

179 

-3 


®This analysis was performed using Lsndsst data from frames 1822-15315 and 
15322. 

■/) - ' ; 
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Total Phosphorus Estimation 

The multiple regression analysis yielded the model: 

log (total phosphorus) = 6.889 ■ 0,060 RB4 
■6 

The model explains about 67 percent of the variance about the mean, 
table 18 . The observed and predicted total phosphorus values, along 
with their residuals, are found in table 19. 


Table 18 

ANALYSIS OF VARIANCE - TOTAL PHOSPHORUS^ 




Analysis 

of Variance 


Source 

DF 

Sum of Squares 

Mean Square 

F 

Regression 

1 

3.239 

3,239 

26.86 

Residual 

13 

1.568 

0.121 

(R2 = 0.674) 

Total 

14 

4.807 


®This analysis 
and 15322. 

was performed using Landsat data 

from frames 

1822-15315 



Table 19 

TOTAL PHOSPHORUS RESIDUALS® 


Reservoir 

Name 

Identification 

Number 

Observed 

Total 

Phosphorus 

(mg/1) 

Predicted 

Total 

Phosphorus 

(mg/1) 

Residual 

Observed-Predicted 

(mg/1) 

Allatoona 

GA 1 

0,026 

0.034 

-0.008 

Blue Ridge 

GA 4 

0.012 

0.015 

-0.003 

Burton 

GA 6 

0.008 

0.010 

-0.002 

Chatuge 

GA 7 

0.016 

0.025 

-0.009 

Cherokee 

TN 9 

0.068 

0.044 

0.024 

Chickamauga 

TN 10 

0 . 032 

0.042 

-0.010 

Douglas 

TN 13 

0.038 

0.042 

-0.004 

Fontana 

NC 14 

0.022 

0.017 

0.005 

Fort Loudoun 

TN 15 

0.060 

0.043 

0.017 

Hiwassee 

NG 18 

0.021 

0.021 

0.000 

Junaluska 

NG 19 

0.035 

0.023 

• 0.012 

Nottely 

GA 22 

0.018 

0.022 

-0.004 

Santeetlah 

NC 27 

0.013 

0.016 

-0.003 

Sidney Lanier 

GA 28 

0.026 

0.013 

0.013 

Watts Bar 

TN 32 

0.033 

0.0.38 

-0.005 


^'ihis analysis was performed using Landsat data from frames 1822-15315 and 
15322. 
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Trophic State Index Estimatioii 

The best reg<rc88ion model foe the pi^ediction of trophic state iadex 
(PCI) values of the 15 NES siiinpled reservoirs using tandsat data from 
frames 1822-15315 and 15322 is ! 

PCI s 4,344 - 0,53‘j MB4 + 0.633 MBS - 2.511 VB7 

The model explains about 88 percent of the variance about the mean, table 20» 


Table 20 

ANALYSIS OF VARIANCE OF REGRESSION MODEL FOR THE 
PREDICTION OF THE TROPHIC STATUS OF 15 RESERVOIRS 
FOUND IN LANDSAT FRAMES 1822-15315 AND 15322 


Analysis of Variance 

Source DF Sum of Squares Mean Square F 


Regression 3 5.826 1.942 27.23 

Residual 11 0.785 0.071 

Total 14 6.611 (r 2 a 0.8813) 


Keep in mind in examining the above model that the trophic state 
index (PCI), as well as values for Secchi disc depth, conductivity, and 
total phosphorus, was developed using mean values of the ground truth 
measurements taken on three occasions during 1973 while the Landsat 
data were coUected within a few seconds on October 23, 1974, 

Figure 14 presents a plot of the observed versus the predicted PCI 

values, Those reservoirs plotted to the upper left of the diagonal have 

been predicted to be more eutcophic than their PCI values indicate. 

Reservoirs plotted to the lower right of the diagonal are estimated to 

be in better condition and less eutrophic than their PCI values suggest. 

When the model was used to predict trophic state index values of 
reservoirs using Landsat data from another date it produced poor results. 
However, it is not unreasonable to expect that over a period of time a 
reservoir would present average reflectance values, which would be more 
representative of its trophic status. This is one area in which it is 
recommended that further investigation be performed, The use of reflec- 
tance values averaged over several times of the year and corrected for 
atmospheric differences and solar angle effects, may result in models 
which could be applied to any reservoir on any given date. 
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GENERATION OF THEMATIC REPRESENTATIOK 


In order to visually represent the relative trophic state of the 
reservoirs I a pixel-by*pixel classification of the Landsat imagery was 
performed, A table- look-up elliptical classifier with a minimum distance 
to mean option was chosen as the algorithm to perform the classification. 

Subsets of landsat scenes 1822-1S31S and 1822-15322 were extracted 
from the computer tapes which included data from each reservoir to be 
classified. In all, 28 subsets were extracted and copied onto another 
tape. A statistics file was generated from the data used in the 
regression models used to predict the trophic classes. From this 
statistics file a table was generated to be used by the classifier. 
Fourteen different classes were represented in the table. 

Each of the subsets of the landsat data were separately classified 
and files with the results stored on tape. To more easily produce the 
thematic representation of the reservoirsi the classified data from each 
scene was combined into two large files. Color separations of each file 
were produced on an Optronics P-1700 photowrite unit. Some classes were 
combined to produce a final thematic representation is six colors (figures 
15 and 16). 

As can be seen in figures 15 and 16, there are variations in the 
color, representing trophic state, within some reservoirs. These 
variations do exist, but may be partially due to nonuniformity in the 
calibration among detectors on the landsat satellite, The thematic 
representations are not as good a representation of the trophic state 
of the reservoirs as the output from the regression models. 



Figure 15. Color coded thematic representation of the trophic status 
of selected reservoirs. 
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Figure 16. Color coded thematic representation of the trophic status 
of selected reservoirs. 

(Landsat scene 1822 • 15315, 23 October. 1974) 



CONCLUSIONS 


This study has shown that there is a definite correlation between 
Landsat reflectance values and reservoir water quality. Correlations 
between MSS imagery and trophic indicators of reservoirs were often 
found to exceed 0.95 in averal Landsat scenes. This correlation has 
been demonstrated wlthoiic the best experimental conditions (no Landsat 
data within one year of ground samples), further indicating the validity 
of this technique. 

It has been shown that a generalized model can be defined to predict 
whole lake trophic status (PCI value) trom Landsat MSS reflectance values. 
Models developed for each of the four Landsat scenes to predict trophic 
status were all significant at the 0.05 level and in all cases had 
values equal to or g’*eater than 0.85. Models were also defined to predict 
specific water quality parameters of Secchi disc depth, conductivity, and 
phosphorus content. Finally, color-coded thematic maps were produced 
showing some of the spatial water quality variance within and between 
reservoirs. 


RECOHMEKDATION S 


I^^ls ifcconunendod this study bu continued so that the following 
concotns' can be further evaluated. 

1. Techniques lor correcting Landsat data for the effects of 
atmospheric absorption » scattering, and sun angle changes 
should be evaluated. If adequate corrections can be made, 
models may be developed for more than one date or scene, or 
the some model may be applicable to other dates or scenes. 

Data has been received to make this evaluation for seven 
scenes using the ground samples and models presented horein. 

2. It is anticipated that much of the variance in these models 
could be eliminated if ground samples are taken concurrently 
with Landsat overpasses. This would require a well coordinated 
sampling program but it should be tried on a demonstration 
basis. 

3. Multiple observations during the growing season have i«en 
valuable in other studies (Scarpace ct al., 1978) to pw-dict 
lake trophic status and type of water quality problem, it is 
recommended that a predictive model which uses multiple data 
of Landsat imagery during one growing se.ison be evaluated, 

4. The primary objective of this study has been to predict whole 
reservoir trophic status. Further studies should also be made 
to name and present the water quality changes within particular 
reservoir areas. 

5. In future studies the design of a ground sampling program should 
strive to produce a more normal di.stribution. This study was 
biased toward eutrophic reservoirs, 

6. The rolatiouship between landcover and water quality should be 
investigated using Landsat data. 

7. Reservoir data extractiou teclmiques that include geometric 
overlays based or actually known reservoir coufigurations should 
be investigated. 
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APPENDIX A 


Correlation Coefficients and Regression 
Models for Landsat Scenes: 

1084-15431 (Oct. 15, 1972), Section 1 
1948-15264 (Feb. 26, 1975), Section 2 
2224-15303 (Sep. 3, 1975), Section 3 


Se^tionJ.. Scoac' 1084-15431 (October 15, 1972) 

Seven reseevoirs were extracted from the frame: 
Boone (TN5) 

Cherokee (TN9) 

Douglas (TN13) 

Fontana (NC14) 

Fort Loudoun (TN15) 

Santeetiah (NG27) 

Watts Bar (TN32) 


Table l.l. Gorrelation between Ground Truth Water 
Quality Data and Landsat Data for Seven 
Reservoirs in Frame 1084-15431 
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PCI 

S^cchl 

Conductivity 

Phosphorus 

Band 4 (MB4) 

0.558 

-0.812 

0.705 

0.458 

Band 5 (MB5) 

0.519 

-0.844 

0.631 

0.441 

Band 6 (MB6) 

0.263 

-0.602 

0.354 

0.178 

Band 7 (MB7) 

-0.405 

0.229 

-0.409 

-0.447 

Band 4/Band 5 (RB4) 

-0,441 

0.806 

-0.519 

-0.389 

Band 5/Band 6 (RB5) 

0.807 

-0.960 

0.908 

0.767 

Band 6/Band 7 (RB6) 

0.813 

-0.986 

0.896 

0.737 

Variance Band 4 (VB4) 

0.030 

-0.435 

0.229 

-0.064 

Variance Band 5 (VB5) 

0.044 

-0.454 

0.237 

-0.048 

Variance Band 6 (VB6) 

0.034 

-0.423 

0.222 

-0.051 

Variance Band 7 (VB7) 

-0.156 

0.238 

-0.196 

-0.102 

Variance Band 4/Baiul 5 





(VRB4) 

0.188 

-0.356 

0.181 

0.280 

Variance Band 5/Band 6 





(VRB5) 

-0.843 

0.897 

-0.898 

-0,804 

Vacianco Band 6/Band 7 





(VRB6) 

-0.522 

0.764 

-0.630 

-0.425 


The best regression model for the prediction of trophic state (PCI) 
isi 


PCI s 4.420 - 0.286MIJ6 - 0.012VRB5 


Table 1.2, Analysis of Variance of Regression Model 
For the Prediction of the T^rophic Status of Seven 
Re.servoirs Found in Lands^'t Frame 1084-15431 


Source 

OF 


Analysis of Variance 


Sum of Squares 

Mean Square 

’ "F " 

Regression 

2 

2.898 

1.449 

12.44 

Residual 

4 

0.466 

0.116 


Total 

5 

3.363 


(R^^ s 0.861) 



Table 1.3. PCI 
Found in ] 

Residuals of 
Dandsat Frame 

Seven Reservoirs 
1084-15431 

Reservoir 

Identification 



Residual 

Name 

Number 

Observed 

Predicted 

Observed -Fred i cted 

Boone 

TN5 

0.441 

0.167 

0.274 

Cherokee 

TN9 

0.714 

0.558 

0.156 

Douglas 

TN13 

0.075 

0.507 

-0.432 

Fontana 

NC14 

-0.960 

-0.600 

-0.360 

Fort Loudoun 

TN15 

0.530 

0.362 

0.168 

Santeetlah 

NC27 

-1.216 

-1.353 

-0.137 

Watts Bar 

TN32 

-0.036 

-0.093 

0.057 


I 



The best regtrcssion model for the prediction of Secchi disc depth let 


log^ (Secchi) s 2S.8S4 - 0.094RB6 


Table 1.4. Analysis of Variance of Regression Model for the 
Prediction of Secchi Disc Depth of Seven 
Reservoirs Found in Landsat Frame 1084*15431 


Source 

DF 

Analysis of Variance 


Sum of Squares 

Mean Square 

F 

Regression 

1 

1.487 

1.487 

180.50 

Residual 

5 

0.041 

0.008 


Total 

6 

1.529 


itr s 0.973 


Table 1.5. Secchi 
Reservoirs Found 

Disc Depth Residuals of Seven 
in landsat Frame 1084-15431 


Reservoir 

Name 

Identification 

Number 

Observed 

Residual 

Predicted Observed-Predicted 

h'oone 

TN5 

57.9 

67.2 

- 9.3 

Cherokee 

TN9 

51.3 

48,9 

2.4 

Douglas 

TN13 

57.2 

54.4 

2.8 

Fontana 

NC14 

107.4 

112.4 

- 5.0 

Fort Loudoun 

TN15 

34.4 

35.8 

^ - 1,4 

Santeetlah 

NC27 

133.6 

122.1 

11.5 

Watts Bar 

TN32 

37.8 

36.0 

1.8 


The best tregi'esslon model tor the prediction of condnctivity is; 


log (conductivity) » -11.119 + 0.092RB5 
c 


Table 1.6. Analysis of Variance of Regression Model 
For the Prediction of Conductivity of Seven 
Reservoirs Found in handsat Frame 1084-15A31 


Source 


Analysis of Variance 
Sum of Squares Mean Square 


Regression 

Residual 

Total 


6.019 

1.280 

7.300 


Table 1,7. Cenductivity Residuals of Seven 
Reservoirs Found in Landsat Frame 1084-15431 


Reservoir 

Name 


Identification 

Number 


Observed Predicted 


Residual 

Observed-Predicted 


Boone 

Cherokee 

Douglas 

Fontana 

Fort Loudoun 

Santeetlah 

Watts Bar 



6 - 


Tlu' bcj$t. regression model for the prediction of total phosphorus is; 


logg (Total Phosphorus) s 3.490 - 0.188MB4 - 0.014VRB5 


Table 1.8. Analysis of Variance of Regression Model for the 
Prediction of Total Phosphorus of Seven Reservoirs 
Found in Landsat Frame 1084-15431 


Source 

DF 

Analysis of Variance 

Sum of Squares 

Mean Squara F 

Regression 

2 

2.022 

l.Oll 22.84 

Residual 

4 

0.177 

0.044 , 

Total 

6 

2.200 

(R'^ = 0.919 

Table 1.9. Total Phosphorus Residuals of Seven Reservoirs 


Found in landsat Frame 

1084-15431 


Reservoir 

Nome 

Identification 

Number 

Observed 

Predicted 

Residual 

Observed-Predicted 

Boone 

TN5 

0.059 

0.053 

0.006 

Cherokee 

TN9 

0.068 

0,059 

0.009 

Douglas 

TN13 

0.038 

0.046 

-0.008 

Fontana 

NG14 

0.022 

0.026 

-0.004 

Fort Loudoun 

TN15 

0.060 

0.067 

-0.007 

Santeetlah 

NC27 

0,013 

0.013 

0.000 

Watts Bar 

TN32 

0.033 

0.026 

0.007 
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Section 2 . Scene 1948*15264 (February 26, 1975) 

Eiglit reservoirs were extracted trom the frame: 
Cherokee (TN9) 

Douglas (TN13) 

Fontana (NC14) 

Fort Loudoun (TN15) 

Junaluska (NC19) 

Santeetlah (NC27) 

WaterviXle (NG31) 

Watts Bar (TN32) 


Table 2.1. Correlation between Ground Truth Water 
Quality Data and Landsat Data for Eight 
Reservoirs in Frame 1948-1526/| 


lOBg 



PCI 

Secchi 

Cdhductivity 

Phosphurus 

Band 4 (MB4) 

-0.070 

0.013 

0.026 

-0.08? 

-0.0|^ 

Band 5 (MBS) 

0.017 

-0.068 

0.108 

Band 6 (MB6) 

0.298 

-0.339 

0.386 

0.241 

Band 7 (MB7) 

0.607 

-0.572 

0,644 

0/616 

Band 4/Band 5 (RB4) 

-0.098 

0.187 

-0,141 

-0.013 

Band S/Band 6 (RB5) 

-0.532 

0,449 

-0.494 

-0.541 

Band 6/Band 7 (RB6) 

-0.357 

0,276 

-0.342 

-0.423 

Variance Band 4 (VB4) 

-0.549 

0.637 

-0.521 

-0.394 

Variance Band 5 (VB5) 

-0,526 

0,563 

-0.499 

-0.368 

Variance Band 6 (VB6) 

-0.485 

0.547 

-0.414 

-0.333 

Variance Band 7 (VB7) 
Viriance Band 4/Band 5 

-0.671 

0.761 

-0.746 

-0.570 

(VRB4) 

Variance Band 5/Band 6 

-0.556 

0.580 

-0.567 

-0.393 

(VRB5) 

Variance Band 6/Band 7 

-0.516 

0.595 

-0,657 

-0.510 

(VRB6) 

-0.671 

0.580, 

-0.591 

-0.434 
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The best reiressioh model for the prediction of trophic ttste (PCI) ie: 


PCI s 5.859 • 0.186HBA - 0.008VRB5 


Table 2.2. Analysis of Variance of Regression Model for the 
Prediction of the Trophic Status of Eight Reservoirs 
Found in Landsat Frame 1946-15264 


Source 

DF 

Analysis of Variance 


Sum of Squares 

Mean Square 

F 

Regression 

2 

3.716 

1.858 

14.00 

Residual 

5 

0.663 

0.133 


Total 

7 

4.379 


(R^ s 0.849) 


Table 2.3, PCI Residuals of Eight Reservoirs 
Found in landsat Frame 1948-15264 

Reservoir 

Identification 



Residual 

Name 

Number 

Observed 

Predicted 

Observed-Predicted 

Cherokee 

TN9 

0.714 

0.292 


0.422 

Douglas 

TN13 

0.075 

0.086 


-O.Oll 

Fontana 

NC14 

-0.960 

-1.145 


0.185 

Fort Loudoun 

TN15 

0.530 

0.414 


0.116 

Junaluska 

NC19 

0.258 

0.041 


0.217 

Santeetlah 

NC27 

-1.216 

-0.835 


'0.381 

Waterville 

NC31 

1.067 

1.122 


•0.055 

Watts Bar 

TN32 

-0.036 

0.456 


-0.492 


The beet tegcession model for the prediction of Secchi disc depth is: 


logg (Secchi) s -1.844 + 0.ia6MB4 + 0.12SNB6 ♦ 0.011VRB6 


Table 2.4. Analysis of Variance of Regression Model for the 
Prediction of Secchi Disc Depth of Eight Reservoirs 
Found in tandsat Frame 1948-15264 


Analys i s of Variance 

SujD of Squares Mean Square 


1.896 

0.140 

2.036 


32 

11 

'35 

/n2 


Table 2.5. Secchi Disc Depth Residuals of Eight 
Reservoirs Found in Landsat Frame 1948-15264 


Identification 

Number 


Observed Predicted 


Residual 

Observed-Predicted 


51.3 
57.2 

107,4 

34.4 

38.0 

133.6 

31.7 

37.8 


41.6 
51.5 
125.0 
41 .8 
37.2 
117.3 
.33.4 
40.4 


9.7 

5.7 
-17.6 

- 7.4 

0.8 

16.3 

- ti-7 

- 2.6 
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The bcKt regr«!i»ion model for the prediction of conductivity 1*; 


log^, (Conductivity) « 13.808 - 0.0A1RB5 - 0.008VRB5 


table 2.6. Analyaia of Variance of Regression Model for the 
Prediction of Conductivity of Eight Reservoirs 
Found in Landsat Frame 1948-15264 


Source 

DF 

Analysis of Variance 


Sum of Squares 

Mean Square 

F 

Regression 

2 

9.584 

4.792 

40.12 

Residual 

5 

0.597 

0.119 


Total 

7 

10.181 


(R* a 0.941) 


Table 

2.7. Conductivity Residuals 
Found in Landsat Frame 

of Eight Reservoirs 
1948-15264 

Reservoir 

Identification 



Residual 

Name 

Number 

Observed 

predicted 

Observed-Predicted 

Cherokee 

TN9 

260 

149 

111 

Douglas 

TN13 

184 

287 

-103 

Fontana 

NC14 

27 

23 

4 

Fort Loudoun 

TN15 

210 

211 

- 1 

Junaluska 

NC19 

106 

114 

- 8 

Santeetlah 

NC27 

17 

22 

- 5 

Waterville 

NC31 

657 

616 

41 

Watts Bar 

TN32 

176 

178 

- 2 


u 


The best regression model for the prediction of totel phosphorus is: 


logg (Total Phosphorus) « 3.521 - 0.238MB4 t 0.091VB4 - 0.009VRB5 


Table 2,8. Analysis of Variance of Rogrcssion Model for the 
Prediction of Total Phosphorus of Eight Reservoirs 
Found in fandsat Frame 1948>15264 


Source 

DF 

Analysis of Variance 


Sum of Sfjuaros 

Mean Square 

j? 

Regression 

3 

2.790 

0.930 

11.75 

Residual 

4 

0.317 

0.079 


Total 

7 

3.107 


(R'‘ s 0,898) 


Table 2,9, Total Phosphorus Residu,als of .iight Reservoirs 
Found in lisodsat Frame 1948*l526u 


Reservoir 

Name 


Identification Residual 

Number Observed Predicted Observed"Predlcted 


Cherokee 

TN9 

0.068 

0.055 

0,013 

Douglas 

TN13 

0.038 

0.038 

0.000 

Fontana 

NC14 

0,022 

0.023 

-0.001 

Fort Loudoun 

TN15 

0.060 

0.048 

0.012 

Junaluska 

NCI 9 

0.035 

0.031 

0.004 

Santcetlah 

NC27 

0.013 

0.014 

-0.001 

Waterville 

NC31 

0.108 

0,313 

-0.005 

Watts Bar 

TN32 

0.033 

0.051 

-0,018 
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Section 3 t Scene 2224" 15303 (September 3| 1975) 

Seven reservoir! were extracted from the frame: 
Boone (TN5) 

Cherokee (TN9) 

Douglas (TN13) 

Fontana (NC14) 

Fort Loudoun (TN15) 

Santeetlah (NC27) 

Watts Bar (TN32) 


Table 3tl. Correlation between Ground Truth Water 
Quality Data and Landsat Data for Seven 
Reservoirs in Frame 2224-15303 



PCI 

Secchi 

logg 

Conductivity 

Phosphorus 

Band 4 (MB4) 

0.328 

-0.619 

0.512 

0.215 

Band 5 (MB5) 

0.88B 

-0.816 

0.895 

0.826 

Band 6 (HB6) 

0,596 

-0.793 

0.746 

0.492 

Band 7 (MB7) 

0.631 

-0.4' 4 

0.582 

0.627 

Band 4/Band 5 (RB4) 

-0.788 

0.941 

-0.891 

-0.728 

Band 5/Band 6 (KBS) 

0.349 

-0.557 

0.520 

0.230 

Band B/Band 7 (RB6) 

-0.160 

-0.080 

-0.048 

-0.220 

Variance Band 4 (VB4) 

0.018 

-0.420 

0.208 

-0.072 

Variance Band 5 (VB5) 

0,018 

-04/. 

0.209 

-0.073 

Variance Band 6 iVfi6) 

0.015 

-0.418 

0.206 

-0.075 

Variance Band 7 (VB7) 

-0.896 

0.669 

-0.872 

•0.872 

Variance Band 4/Band 5 





(VRB4) 

-0.219 

-0.223 

0.006 

-0,310 

Variance Band 5/B<f il 6 





(VRB5) 

-0.646 

0.256 

-0.548 

-0.672 

Variance Band B/Band 7 





(VRB6) 

-0.674 

0,640 

-0.749 

-0.585 


The best regression model for the prediction of trdphic state (PCI) is: 

PCI s 6.127 + 0.006VRB5 - 11.127VB7 , 

/ 

/ 

Table 3.2. Analysis of Variance of Regression Model for the 
Prediction of the Trophic Status of Seven Reservoilrs 



Found 

in Landsat Frame 2224-’15303 { 

1 





Analysis of Variance 


Source 

DF 

Sum of Squares 

Mean Square 

F 

Regression 

2 

2.972 

1.486 

I' 

15,20 

Residual 

4 

0.391 

0.098 

9 

Total 

6 

3.363 


(R^ = 0.884) 



Table 3.3. PCI Residuals of 
Found in Landsat Frame 

Seven Reservoirs '( 

2224-15303 

Reservoir 

Identification 



Residual 

Name 

Number 

Observed 

Predicted 

Observed-Predicted 

Boone 

TN5 

0.441 

0.254 

0.187 

Cherokee 

TN9 

0.714 

0.425 

0.289 

Douglas 

TN13 

0.075 

0.562 

-0.487 

Fontana 

NC14 

-0.960 

-1,045 

0.085 

Fort Loudoun 

TN15 

0.530 

0.453 

0.077 

Santeetlah 

NC27 

-1.216 

-1.069 

-0.147 

Watts Bar 

TN32 

-0.036 

-0.031 

-0.005 





1 





The best regression model for the prediction of Sccchi disc depth is: 


log^ (Sccchi) » «8.959 + 0,083RBA 


Table 3.4. Aluilysis of Variance of Regression Model for the 
Reduction of Sccchi Disc Depth of Seven Rescrvoi .’s 
Found in Landsat Frame 2224-15303 


Source 

1)F 

Analysis of Variance 


Sum of Squares 

Mean Square 

F 

Regre.ssion 

1 

1.354 

1.354 

38.80 

Residual 

5 

0.175 

0.035 


Total 

6 

1.529 


(R*^ =! 0.886) 


Table 3.5. Sccchi Dise Dopth Residuals of Seven Reservoirs 
Found in Landsat Frame 2224-15303 

n i 


Reservoir 

Name 

Identification 

Number 

Observed 

Predicted 

Residual 

Observed-Predicted 

Boone 

TN5 

57.9 

68.5 

-10.6 

Cheriikee 

TN9 

51.3 

55.6 

-4.3 

Douglas 

TN13 

57.2 

41.0 

16.2 

Fontana 

NC14 

107.4 

99.9 

7.5 

Fort Loudoun 

TN15 

34.4 

36.5 

- 2.1 

Santeetlah 

NC27 

133.6 

128.2 

5.4 

Watts Bar 

TN32 

37.8 

43.4 

- 5.6 



''4 

i 










15 


The best regression model for the prediction of conductivity is; 


log (Conductivity) a 12.002 - 0.00j,vp'v - 11.561VK7 


Table 3.6. Analysi,s of Variance of Regression Model for the 
Prediction of Conductivity of Seven Reservoirs 
Found in Landsat Frame 2224*^15303 


Source 

DF 

Analysis of Variance 

Sum of Squares 

Mean Square F 

Regression 

1 

7.028 

3.514 51.84 

Residual 

4 

0.271 

0.068 5 

Total 

6 

7.299 

(R*^ s 0.963) 


Table 3,7. Conductivity Residuals of Seven Resorvoir.s 
Found in bandsat Frame 2224-15303 


Reservoir 

Name 

Identification 

Number 

Observed 

Predicteil 

Residiu 

Observed-Prc 

Boone 

TN5 

174 

143 

31 

Cherokee 

TN9 

260 

241 

19 

Douglas 

TNI 3 

184 

268 

-84 

Fontana 

NC14 

27 

24 

3 

Fort Loudoun 

TN15 

210 

176 

34 

Santeetlah 

NC27 

17 

21 

- 4 

Watts Bar 

TN32 

l7o 

176 

0 
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The best regressiaii moriel for the prediction of totnl phosphorus isi 


loBj, (Total Phosphorus) « 0.324 - 5.S69VB7 


Table 3.8. Analysis of Variance of Regcesston Model for the 
Prediction of Total Pliosphorus of Seven Reservoirs 
Found in landsat Frame 2224-15303 


Analys i s of Variance 


Source 

DF Sum 

of Squares 

Mean Square 

" ' F 

Regression 

1 

1.671 

1.671 

15.84 

Residual 

5 

0.528 

0.106 


Total 

6 

2,199 


(R'^ = 0,760) 

Table 

3.9. Total Phosphorus Residuals of Seven Reservoirs 


Found in Lnndsat Frame 

2224-15303 


Reservoir 

Identification 



Residu.ll 

Name 

Number 

0b.secved 

Predicted Observed-Predicted 

Boone 

TN5 

0.059 

0.047 

0.012 

Cherokee 

TN9 

0.068 

0,060 

0.008 

Douglas 

TN13 

0.038 

0 . 063 

-.0.025 

Fontiina 

NCI 4 

0 . 022 

0.020 

0.002 

Port Loudoun 

TNI 5 

0.060 

0.051 

0.009 

Santcctlah 

NC27 

0.013 

0.018 

-0.005 

Watts Bar 

TN32 

0.033 

O’. 027 

0.006 


C^<c' 


